function [DWL_pertax,DWL_perinc]=show_DWL(p_interest_univ,y_interest_univ,p_interest,y_interest,tau_vec,K1,K2,K3,p_lb,p_ub,q_lb,q_ub,y_lb,y_ub,mu1,mu2,mu3,sigma2,sigma3,fdist1,fdist2,fdist3,alpha1,alpha2,alpha3,est,q_hat_est,show_txt,addtxt,fname,flag_median)

	p_interest_univ_log=p_interest_univ;
	y_interest_univ_log=y_interest_univ;
	y_interest_log=y_interest;

	DWL=NaN*ones(size(tau_vec,1),3); 
	DWL_pertax=DWL;
	DWL_perinc=DWL;
	DWL_tri_pertax=DWL;
	DWL_tri_perinc=DWL;

	tt_max=size(tau_vec,1);
	tt_min=1;

	if flag_median==1
		tt_min=2;
		tt_max=2;
	end

	for tt=tt_max:(-1):tt_min
	
		q_interest=q_hat_est(:,tt); 
		
		[~, ~, auxBB_deriv2_interest, auxBB_deriv3_interest] = construct_chebyshev_berkson_Ponly_twodim_inc(p_interest,q_interest,y_interest,K1,K2,K3,p_lb,q_lb,y_lb,p_ub,q_ub,y_ub,NaN,NaN,NaN,NaN,NaN,NaN,mu1,mu2,mu3,0*p_interest,sigma2,sigma3,fdist1,fdist2,fdist3,alpha1,alpha2,alpha3);

		aux_Ginv_y_interest_est=auxBB_deriv3_interest*est.theta_hat; 
		aux_Ginv_q_interest_est=auxBB_deriv2_interest*est.theta_hat; 
		aux_dg_dy_interest_est=-aux_Ginv_y_interest_est./aux_Ginv_q_interest_est;
		aux_g_level_est=exp(q_hat_est(:,tt));
		aux_g2_level_est=aux_dg_dy_interest_est.*aux_g_level_est./exp(y_interest);

		for yy=1:3

			[DWL(tt,yy), DWL_pertax(tt,yy), DWL_perinc(tt,yy), ~, DWL_tri_pertax(tt,yy), DWL_tri_perinc(tt,yy)]=compute_DWL(p_interest_univ_log,y_interest_univ_log,y_interest_log,aux_g_level_est,aux_g2_level_est,yy);

			if show_txt==1
				fileID=fopen(fname,'a');
				fprintf(fileID,'%5.3f; %6.0f; %6.3f; %4.2f; %s\n',tau_vec(tt),exp(y_interest_univ_log(yy)),DWL_pertax(tt,yy),10^4*DWL_perinc(tt,yy),addtxt);
				fclose(fileID);
			end

		end

	end
	
	%

